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ABSTRACT 

The aberrated radiation pressure at the inner edge of the accretion disk around an astrophysical black 
hole imparts a relative azimuthal velocity on the electrons with respect to the ions which gives rise 
to a ring electric current that generates large scale poloidal magnetic field loops. This is the Cosmic 
Battery established by Contopoulos and Kazanas in 1998. In the present work we perform realistic 
numerical simulations of this important astrophysical mechanism in advection-dominated accretion 
flows-ADAF. We confirm the original prediction that the inner parts of the loops are continuously 
advected toward the central black hole and contribute to the growth of the large scale magnetic 
field, whereas the outer parts of the loops are continuously diffusing outward through the turbulent 
accretion flow. This process of inward advection of the axial field and outward diffusion of the return 
field proceeds all the way to equipartition, thus generating astrophysically significant magnetic fields 
on astrophysically relevant timescales. We confirm that there exists a critical value of the magnetic 
Prandtl number between unity and 10 in the outer disk above which the Cosmic Battery mechanism 
is suppressed. 
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1. INTRODUCTION 

The origin of astrophysical magnetic fields remains an 
open issue of modern astrophysics. Ludwig Biermann, in 
his famous paper, proposed a mechanism for the genera¬ 
tion of large-scale electric currents based on the thermo¬ 
electric effect (Biermann 1950). This mechanism leads 
to magnetic fields that are usually quite weak initially 
and require subsequent dynamo amplification to reach 
astrophysically relevant magnitudes. Some years ago we 
proposed a mechanism alternative to that of the Bier¬ 
mann Battery (Contopoulos & Kazanas 1998, hereafter 
CK, Contopoulos et al. 2006), the so called Cosmic Bat¬ 
tery. We posited that the aberration of the radiation 
field experienced by the electrons of an accretion flow 
around an astrophysical source (a generalization of the 
well known Poynting-Robertson effect; Poynting 1903, 
Robertson 1937, Lamb & Miller 1995, Bini et al. 2009, 
Koutsantoniou & Contopoulos 2014) generates toroidal 
electric currents sufficiently large to support poloidal 
magnetic fields that under certain astrophysically plausi¬ 
ble circumstances will grow to equipartition values. The 
aberrated radiation force implies a direct coupling be¬ 
tween accretion flow vorticity, disk plasma diffusivity and 
the large scale ordered magnetic field which may have di¬ 
rectly observable implications (Contopoulos et al. 2009). 

The magnetic field loops generated by the action of the 
Cosmic Battery around the inner edge of the accretion 
disk (where radiation pressure and aberration effects are 
maximal) are anchored on different radii of the inner ac¬ 
cretion disk, thus they become twisted in the azimunthal 
direction by the differential rotation of the flow. As their 
twisting relaxes in the vertical direction, the loops open 
up and separate into an inner component and an outer 
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component (the return field) that threads the disk fur¬ 
ther away from the axis. The poloidal fields of the two 
components are in opposite directions, the inner one is 
parallel and the outer one is antiparallel to the angular 
velocity vector. On top of that, an overall helical mag¬ 
netic field is produced when the footpoints of the initially 
poloidal magnetic field loops are dragged by the rotat¬ 
ing disk plasma, with the inner toroidal magnetic field 
pointing in the direction opposite to that of the disk ro¬ 
tation in the northern hemisphere of the disk, and along 
the direction of the disk rotation in the southern hemi¬ 
sphere. According to the Cosmic Battery, reversing the 
observer’s hemisphere, or equivalently, the direction of 
disk rotation, reverses the polarity of the axial field but 
leaves the direction of the toroidal field unchanged in 
the observer’s sky. This can be tested in radio VLBI jets, 
where the direction of the toroidal field can be directly in¬ 
ferred from measurements of Faraday rotation gradients 
transverse to the jet axis (Blandford 1993, Reichstein & 
Gabuzda 2012, Murphy & Gabuzda 2013). Observations 
of the toroidal fields revealed from such measurements 
in pc-scale jets from active galactic nuclei show that, in 
most cases, the axial electric current in the VLBI core jet 
flows toward the center (Contopoulos et al. 2009). The 
probability that this asymmetry came about by chance 
was calculated to be less than 1%. This observational re¬ 
sult supports the hypothesis that the universe is seeded 
by cosmic magnetic fields that are generated in active 
galactic nuclei via the mechanism of the Cosmic Battery, 
and are subsequently injected in intergalactic space by 
the jet outflows. 

In the original Cosmic Battery paper (CK) we studied 
the Poynting-Robertson effect of a central isotropic radi¬ 
ation field on the evolution of the magnetic field in a thin 
accretion flow. We followed the evolution of the magnetic 
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field by solving the induction equation in ID under the 
assumption that the generated magnetic field is purely 
vertical. We extended this analysis in axisymmetric con¬ 
figurations (2D; Contopoulos et al. 2006, Christodoulou 
et al. 2008) under the assumption that the accretion disk 
is connected to a steady-state force-free magnetosphere 
above and below it. Our work faced some criticisms. 
The main argument against the Cosmic Battery con¬ 
sists of the misconception that the poloidal magnetic field 
loops advected to the center saturate the magnetic field 
to values well below equipartition (Bisnovatyi-Kogan et 
al. 2002). This criticism does not take into account the 
fact that because the loops are stretched by the differen¬ 
tial rotation, they open up vertically, and only the inner 
part of the loop is advected inward. In most cases of 
astrophysical significance, the outer part of the loop dif¬ 
fuses outward because astrophysical accretion disks are 
turbulent and viscous, thus also diffusive, and the field 
growth does not saturate. 

We decided to extend our previous axisymmetric sim¬ 
ulations to realistic astrophysical scenaria. We consider 
general hot accretion flows around a spinning black hole. 
The main new element of our present simulations is that 
we now follow the time-dependent evolution of the force- 
free magnetospheric field. As we will see, the magne- 
tospheric field is twisted by the differential rotation in 
the disk and the central spinning black hole. The inner 
part is advected toward the black hole horizon, and the 
outer part diffuses through the outer accretion disk. We 
obtained several animations (movies) of the field evolu¬ 
tion which the interested reader may find in the web (e.g. 
http://youtu.be/keF0V8xwons). Our present analysis 
is Newtonian (i.e. non-general relativistic), but we ex¬ 
pect it to be also valid qualitatively when generalized in 
general relativity. 

In § 2 we present a brief overview of various disk mod¬ 
els, and outline the main physical processes that we con¬ 
sider in our integration of the induction equation. In 
§ 3 we discuss how we couple the disk with a force-free 
magnetosphere above and below it. In § 4 we discus the 
fundamental importance of the magnetic Prandtl num¬ 
ber with two characteristic examples, one that exhibits 
continuous magnetic field growth, and one that exhibits 
saturation of the mechanism. Finally, in § 5 we discuss 
various astrophysical examples of ADAF disks and reach 
our conclusions. 

2. ADAF DISKS AND THE INDUCTION EQUATION 

Black hole accretion is a physical process that pro¬ 
vides the primary source which powers active galactic 
nuclei (hereafter AGN), black hole X-ray binaries (here¬ 
after XRB) and possibly also gamma-ray bursts (here¬ 
after GRB). Black hole accretion flows are divided into 
two categories, cold and hot. Cold accretion flows are 
described by the standard thin disk model (Shakura & 
Sunyaev 1973, Novikov & Thorne 1973, Lynden-Bell & 
Pringle 1974). In this case the disk is geometrically thin 
but also optically thick and quite cold. The thin disk 
model is applicable when the mass accretion rate is some¬ 
what below the Eddington rate. In this paper we study 
the cosmic battery in the general case of a hot accretion 
flow around a rotating black hole. 

The first model for hot accretion was obtained by 
Shapiro et al. (1976). In this case the plasma is more 


tenuous and hotter than the plasma in a thin disk, and 
is also optically thin. The main characteristic of this 
model is the introduction of a two-temperature accret¬ 
ing plasma where the ions have higher temperature than 
the electrons. This model succeeded for the first time in 
explaining the hard X-ray emission seen in black hole 
sources. Unfortunately, it is thermally unstable, and 
therefore it cannot be a good representation for astro- 
physical accretion disks. Ichimaru (1977) was the first 
who studied the important role of advection in hot ac¬ 
cretion flows. He pointed out that, in certain regimes, 
the viscously dissipated accretion energy can go prefer¬ 
entially into heating the accretion flow than into radia¬ 
tion (Ichimaru 1977, Rees et al. 1982, Yuan & Narayan 
2014). This is the central element of advection domi¬ 
nated accretion flows (hereafter ADAF). The discovery 
of ADAF solutions (Narayan & Yi 1994, Narayan & Yi 
1995a, Narayan & Yi 1995b, Abramowicz et al. 1995, 
Chen et al. 1995, Yuan & Narayan 2014) gave rise to 
the application of ADAF models in various black hole 
systems. These systems include the supermassive black 
hole in the Galactic center, Sagittarius A, low-luminosity 
AGNs and black hole X-ray binaries (e.g. Narayan et al. 
1998, Abramowicz & Fragile 2013, Narayan & McClin- 
tock 2008, Ho 2008, Quataert 2001, Lasota 1999). 

ADAF models are characterized by the following six 
equations for the radial velocity u r , the azimuthal veloc¬ 
ity Vrj ,, and the isothermal sound speed c s (Narayan & 
Yi 1994), 
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where, Vk = {GM /r ) 1 / 2 is the Keplerian velocity, M is 
the mass of the central compact object/black hole, and 

c s is the speed of sound, e = where 7 is the ratio 

of specific heats. 7 = | in a relativistic plasma. The ad¬ 
vection parameter 0 < / < 1 measures the fraction of the 
viscously generated heat that is advected with the flow. 
When / —> 1, most of the generated heat is advected with 
the flow, whereas when / —> 0 the accretion flow radi¬ 
ates away all of its generated heat. In that case we end 
up with a standard thin Shakura-Sunyaev disk (Shakura 
& Sunyaev 1973) with accretion constant a ~ 0.1 (this 
value corresponds to accretion disks around black holes; 
Komissarov 2006). We will henceforth work in spherical 
coordinates (r,9,(f>). Notice that the expression for the 
sound speed is needed in order to determine the scale 
height of the disk h = rc s /vx , or equivalently the polar 
angle of the surface of the disk 
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Eqs. m are non-relativistic. On the other hand, the 
Cosmic Battery is most effective around the inner edge 
of the accretion disk around an astrophysical black hole 
which is believed to coincide with the radius of the in- 
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nermost stable circular orbit (hereafter ISCO). We do 
not plan to perform a full general relativistic numerical 
simulation in that region. Instead, we extend our New¬ 
tonian formulation down to the black hole horizon. We 
believe that, the exact details of the accretion flow in 
that region are not particularly important. As we will 
see below, what is most important is the fact that inside 
the ISCO, the flow transitions from being mostly Keple- 
rian, thus also turbulent, viscous, and diffusive, to being 
mostly freely falling, thus also laminar, non-viscous, and 
non-diffusive. 

Our aim in this paper is to investigate the action of the 
Cosmic Battery in realistic astrophysical accretion disks. 
We will thus assume the above velocity profile inside the 
disk (eqs. 1) and integrate the induction equation (eq. [3] 
below) in order to follow the evolution of the disk mag¬ 
netic field in the presence of the Cosmic Battery radiation 
source. In other words, we will assume that the field is 
passive and does not affect the velocity field in the disk. 

The induction equation in the disk (0disk < 0 < 90°) 
can be written as 

— 7 ^^- = —cV x (^Ecb4> - — x B disk + r/V x Bdisk) 

( 3 ) 

where Bdisk is the magnetic field, and 77 is the magnetic 
diffusivity in the disk. The Cosmic Battery electromo¬ 
tive source term Ecb is equal to the aberrated radiation 
force on the electrons in the azimuthal direction TladU 
divided by the proton charge e. In our original paper, we 
considered a simplified central radiation field with lumi¬ 
nosity L , in which 
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where, <tt is the electron Thompson cross section. The 
correct expression for the aberrated radiation force re¬ 
quires a complex general relativistic calculation of the 
radiation field as felt by the moving electron, taking into 
consideration the plasma optical depth and the fact that 
the source of radiation is the accretion disk itself. This 
calculation has been performed for the first time by Kout- 
santoniou & Contopoulos (2014) for an optically thick 
flow with a thermally emitting surface. The case of an 
optically thin disk is much more complex and is currently 
under investigation. In order to proceed, we will assume 
a slightly more general phenomenological expression 
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which incorporates the main features of a realistic ra¬ 
diation force: the buildup of the radiation field as we 
approach the ISCO, the Poynting-Robertson (radiation 
aberration) effect, and the 1 /r 2 drop-off with distance 
of the radiation pressure in an optically thin medium. 
We also have the option to account for some long op¬ 
tical depth (on the order of nsco) in the disk (second 
expression in top line of eq. [5] above). As we said above, 
a more detailed general relativistic calculation of the ra¬ 
diation force that takes into account the physical optical 
depth of the accretion flow is under way (Koutsantoniou 


& Contopoulos in preparation). 

The magnetic diffusivity 77 can be expressed as 

{ QI ~c 2 . c ^ 

^ lf r - nsco 

( 6 ) 

0 if r < nsco 

(see for example Lovelace et al. 2009). The magnetic dif¬ 
fusivity is therefore a function of the magnetic Prandtl 
number V m which itself is in general a function of dis¬ 
tance in the accretion flow. In the present work we will 
assume that V m is roughly a constant within the spatial 
extent of our numerical simulations, but very quickly be¬ 
comes infinite (i.e. 77 drops very quickly to zero) inside 
nsco as the flow very quickly changes character from 
quasi-Keplerian to almost free fall on the black hole hori¬ 
zon. This transition is related to the fact that as the flow 
rotation decreases, the shear decreases, and the flow be¬ 
comes laminar inside the ISCO because matter is almost 
in free fall there. The laminar flow implies that the tur¬ 
bulence responsible for both the flow viscosity and the 
magnetic diffusivity ceases. 

2 .1. Numerical setup 

We solve our equations in a numerical grid with uni¬ 
form latitudinal spacing 59, linearly increasing radial 
spacing Sr = 3 r56, and a finite radial extent 7 m < r < 
r 0 ut, where r in = ttsco/3, and 9 min <8 < 90°. Although 
our formalism is Newtonian, the inner radial boundary 
ri n of our grid is taken to represent the black hole hori¬ 
zon. An inner polar cutoff of a few degrees 0 m ; n ~ 5° 
is implemented in order to avoid coordinate system sin¬ 
gularities around the axis. Our standard r x 9 grid res¬ 
olution is 87 x 90 respectively, and our numerical grid 
extends radially out to r ou t = 25rigco- We have also run 
selected simulations with double grid resolution, namely 
174 x 180, and found agreement within less than 5%. 
We have thus concluded that our present numerical re¬ 
sults have converged. We have also run simulations on 
a uniform (Sr, 59) grid and obtained essentially the same 
results. 

The magnetic field is generated by the Cosmic Battery 
and is not brought in from large distances, therefore, the 
time integration of eq. ([3]) begins with 

Bdisk(t = 0) = 0 . (7) 

Also, we set B r = = 0 along the equator 9 = 90°, 

and implement free inner and outer boundaries at r- m and 
r 0 ut respectively. What is most important is that the 
magnetic field along the surface of the disk is matched 
to the magnetospheric field, namely 

Bdisk(^disk) = BMs(^disk) ■ (3) 

The magnetospheric field Bms(a 0) for 0 < 0disk is ob¬ 
tained in an independent calculation that we are now 
going to discuss. 

3. THE FORCE-FREE MAGNETOSPHERE 

Whenever the magnetic field generated by the Cosmic 
Battery in the disk reaches the surface it will extend 
into the disk atmosphere, the so-called magnetosphere. 
In order to follow the evolution of the magnetic field in 
the accretion flow, we must also numerically evolve the 
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Tm - 0.4 
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Figure 1 . Poloidal magnetic field lines at times t acc (a), 3t aC c (b), S.6t acc (c), 15.7£ aC c (d), 23.7t acc (e), and 38.4£ acc (f), in the standard 
case of a diffusive accretion disk with Prandtl number = 0.4. Length scales in unit of nsco- Colors represent the value and sign of B$ 
in arbitrary units (blue: negative, orange: positive). The surface of the disk lies at 0 = 60°. Poloidal field loops are continuously generated 
around nsco- The inner footpoints of the field lines are advected together with infalling matter to the central object, and magnetic flux of 
one polarity is accumulated on the central object at a steady rate. The disk is threaded by the return field polarity. The outer footpoints 
are diffusing outwards inside the accretion disk, and because of the differential rotation of their footpoints, the loops open up. A steady 
configuration is established where poloidal loops are continuously generated around the inner edge of the disk and open up to infinity. 
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field in the magnetosphere just above the surface of the 
disk. It is natural to assume that, near the surface of 
the disk, the magnetosphere attains a force-free config¬ 
uration, in analogy to the solar corona and the pulsar 
magnetosphere. The equations that describe the field 
evolution are those of force-free electrodynamics (FFE; 
Gruzinov 1999, Blandford 2002), namely Maxwell’s equa¬ 
tions in the presence of electric charges and currents, 


<9B M s 

dt 


—cV x E M s 


<9E M s 

dt 


= cV x B M s 


-4ttJ 


V ■ Bms =0 


V • EmS =4-7T Pe , 


(9) 


supplemented by the force-free and ideal MHD condi¬ 
tions 


p e EMS + J x Bms = 0 (10) 

Ems • Bms = 0 • (H) 

Here, p e , J are the magnetospheric electric charge and 
current densities respectively. The above formulation is 
obviously non-general relativistic, yet we believe that the 
simulations performed in this study can give an overall 
impression about how the Cosmic Battery mechanism 
works. Notice that recent general relativistic simulations 
of force-free disk magnetospheres (Parfrey et al. 2015) 
show a magnetic held evolution similar to the one we 
obtain here (magnetic held loops becoming twisted and 
opening up due to the differential rotation of the how 
in the disk, current sheet and plasmoid formation, etc.). 
We plan to extend our numerical calculations in general 
relativity in a forthcoming paper. 


3.1. Numerical setup 

The numerical code we employ is the 3D code hrst 
developed in Kalapotharakos & Contopoulos (2009), re¬ 
written in spherical coordinates (r, 9 , <j>) by Contopou¬ 
los (2013). We have lowered the dimensions to 2D by 
assuming axisymmetry. As in the disk interior, we imple¬ 
ment an inner free boundary at r = n n in this region too. 
In practice, we allow for any magnetic flux that is gener¬ 
ated inside the disk and reaches the inner magnetospheric 
boundary (which represents the black hole horizon), to 
be freely distributed along that boundary. Beyond a 
certain distance rpML < r 0u t, we assume an absorbing 
non-rehecting outer boundary in the form of a special 
implementation of an axisymmetric Perfectly Matched 
Layer (PML; discussed in detail in Kalapotharakos & 
Contopoulos 2009), namely 


where 
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( 12 ) 

(13) 


We have taken a 0 = 100, and rpML = 17?tsco- These 
parameters, as well as the exponent in eq. (1131) have been 


chosen empirically. Notice the similarity with the equiv¬ 
alent expression in Cerutti et al. (2014). We have run 
several simulations at various radial resolutions and var¬ 
ious values of the outer boundary r ou t and the results 
are similar (the simulations have converged). For most 
of the results we obtain in the present paper, the outer 
boundary was placed at 25risco- 
As we said in the previous section, the evolution of 
the magnetic fields in the two domains (the disk and 
the magnetosphere) are related. It is important to no¬ 
tice here that there is an important complication that 
enters in the matching between the two domains, and 
this has to do with the abrupt transition between the 
dense turbulent accreting flow and the almost empty ion¬ 
ized (probably also evaporating and outflowing) magne¬ 
tosphere. This is very similar to the well known transi¬ 
tion region in the sun where, within a very thin layer of 
about one hundred kilometers, the temperature rises by 
almost two orders of magnitude to about one million de¬ 
grees, and the matter density drops by about one order 
of magnitude. The solar transition region is the subject 
of ongoing intense theoretical and numerical investiga¬ 
tions. The astrophysical disk transition region is at least 
equally complicated since it is not only the base of the 
disk corona, but also the origin of purported disk winds 
and outflows (e.g. Li 1995, Ferreira & Pelletier 1995). is 
not the aim of the present study to solve in detail for the 
structure of the flow and magnetic field in this region. In 
practice, we define a transition region of two 0-grid zones 
above 0 d isk- I n that region we set 

Bms = B disk . (14) 

This matching allows for magnetic field loops that reach 
the disk surface zones to escape to the magnetosphere. 
We also set the held lines in that region in rotation by 
introducing a magnetospheric poloidal electric Held 

EmS ($disk) =- ~<j> x Bdisk(^disk) • (15) 

c 

Finally, we slightly modify the induction equation that 
we solve for the disk material in that region by introduc¬ 
ing an extra phenomenological term that mimics various 
complex physical effects taking place in the surface layers 
of astrophysical accretion disks such as surface convec¬ 
tion, buoyancy of magnetic Held loops, and the fact that 
the surface layers are highly ionized (thus also fully con¬ 
ducting) due to cosmic ray irradiation, namely 

= "CV X (e C b4> + X B tra ns) , (16) 

where, B tla ns is the disk magnetic field in the. The intro¬ 
duction of this transition layer effectively shields the disk 
interior: a) it prevents any generated magnetospheric 
field from ‘entering’ the disk from above (especially in 
the case of turbulent high p flow conditions), and b) it 
‘pushes outward’ (in the 9 direction) at a very small frac¬ 
tion of the speed of light any magnetic field loop that 
enters it from below. We implemented different fractions 
(1/200, 1/20, etc.) and the results were qualitatively 
similar, transition region. As we said above, we are not 
claiming that we study in detail the disk-magnetosphere 
transition region. This is the reason we opted for a min¬ 
imal two 0-zone layer that only serves the purpose of 
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shielding the disk interior from the disk magnetosphere. 
Notice that Parfrey et al. (2015) ignored the important 
physical significance of a transition layer by directly cou¬ 
pling their magnetosphere to a prescribed distribution of 
the magnetic field on the surface of the accretion disk, 
without solving for the magnetic field distribution in the 
disk interior. 

The inner boundary (which represents the black hole 
horizon) also introduces a rigid body rotation expressed 
through 


EMs(n n , 0 ) = — 


^^(Pin; f^disk) sill(d) 
csin(6 l disk) 


</>xB M s(nn,60 (17) 


for 6 < 0disk- The latter turns out not to be important 
in the determination of the overall magnetospheric field 
topology, and plays a role only if we are interested in the 
extraction of energy from the black hole rotation (Bland- 
ford & Znajek 1977). As in the disk interior, the time 
integration of eqs. (151fl2ll begins with 


B M s(* = 0) = 0 (18) 


everywhere in the magnetosphere. 

4. THE COSMIC BATTERY IN ACTION 


The goal of this paper is to convince the reader that the 
Cosmic Battery plays a fundamental role in the genera¬ 
tion of the large scale magnetic fields that thread astro- 
physical accretion disks and are believed to be respon¬ 
sible for such diverse astrophysical phenomena as jets, 
AGNs, XRBs, and possibly also GRBs. We performed 
numerical integrations that couple the FFE equations in 
the magnetosphere above and below the disk (eqs. lifUTTI) 
with the induction equation that includes the Cosmic 
Battery electromotive force and the flow magnetic dif- 
fusivity in the disk (eq. [3|. As we said, the coupling 
is performed in a thin surface transition region through 
eqs. dump - Our goal is to find out how much magnetic 
flux can be accumulated around the central compact ob¬ 
ject. In particular, we are going to show under what 
conditions the accumulated magnetic flux grows steadily 
to equipartition, and under what conditions it saturates 
to a very low value. We want to make clear that these 
are not MHD flow simulations since the accretion flow is 
independently specified by the ADAF equations (eqs. [T|). 

In this section, we are going to highlight the fundamen¬ 
tal role of the accretion disk magnetic diffusivity in acti¬ 
vating (or under different conditions saturating) the Cos¬ 
mic Battery. This point, already discussed in the original 
CK presentation, has been missed by Bisnovatyi-Kogan 
et al. (2002) in their criticism of the Cosmic Battery. We 
will thus now focus on two particular simulations, one 
with moderate and one with low magnetic diffusivity. 
The first simulation corresponds to a magnetic Prandtl 
number Vm = 0.4. We will call this case the ‘standard’ 
case. In the second simulation we set V m = 100. We 
will call this case the ‘very low rf case. Notice that, in 
order to minimize the number of free parameters in the 
problem, we fix the accretion constant a = 0.1, and the 
advection parameter f = 1 (the latter value corresponds 
to a strongly advection-dominated flow). We will now 
compare the evolution of our two characteristic simula¬ 
tions. 

In the beginning, the Cosmic Battery term (eq. [5J 
generates poloidal magnetic loops around nsco dr both 


cases. The inner footpoints of the field lines are advected 
by disk matter that falls into the central object. Further¬ 
more, as each loop enters the magnetosphere it starts to 
inflate due to the differential rotation between its inner 
and outer footpoints on the surface of the disk. This 
can be seen in Figs, la and 2a which are very similar. As 
the mechanism continues to operate, more and more field 
lines are brought to the central object. The similarities 
end beyond about one accretion time defined as the time 
needed to traverse a distance equal to nsco at the ISCO 
accretion speed, namely 

j. _ nsco 

^acc — / \ 

Ur (nsco) 

_/ 5 + 2e/A nsco 

V 3 a ) UR-(nsco) 

“ L7xltl “ ! (l5Sj) aec - (19) 

What differentiates the two cases is the behavior of the 
outer footpoints beyond that time. 

4.1. The standard case 

In the standard case, the outer footpoint of a poloidal 
field line generated by the Cosmic Battery begins to dif¬ 
fuse outward inside the accretion disk. Because of the 
difference in the rotation between the two footpoints, the 
field lines wind up around the central object. They in¬ 
flate in the magnetosphere above the accretion disk and 
sustain a connection between the central object and the 
disk (Fig. lb). New loops are continuously generated 
around the position of the ISCO (Fig. lb-f). As the outer 
footpoints are continuously moving outward, the loops 
expand more and more. Their inner and outer parts 
never disconnect, and the sign of the magnetospheric B$ 
component remains unchanged (negative). Notice that 
although develops highly variable structure on small 
scales (as can be seen in Figs, le & f), it does not af¬ 
fect the large scale field structure, and the numerical 
integration is stable. The return polarity diffuses out¬ 
ward and the field in the disk remains at very low (sub- 
equipartition) values. Such low magnetic field cannot be 
responsible for a disk MHD wind (e.g. Li 1995, Ferreira 
& Pelletier 1995). This is an interesting point that may 
be investigated more in the future. 

If the central object is rotating, then the field lines of 
one polarity that have already reached it can extract its 
rotational energy (Blandford & Znajek 1977, Nathanail 
& Contopoulos 2014). In this paper, we are interested 
in showing the growth of the magnetic field, and there¬ 
fore, we are not going to discuss this effect any fur¬ 
ther. It is interesting to notice that, for the particular 
ADAF parameters of our present simulation, the inner 
magnetospheric boundary is set in uniform rotation with 
D(r in ) = n disk (rin,0disk) = 0.6c/nsco- This corresponds 
to a magnetospheric light cylinder at cylindrical radius 

Rlc = 777 —t = 1.67nsco ■ ( 20 ) 

S2(n n ) 

As shown in Fig. 3, the magnetic flux brought to the 
central object grows at a steady rate which has no rea¬ 
son to cease before the held reaches equipartition with 
the accretion how (dehned as the limit where the accu- 
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Figure 2. Same as Fig. 1 for the very low 77 case of a non-diffusive optically thin accretion disk with magnetic Prandtl number Vm = 100. 
The inner foot points of the field lines are advected together with infalling matter onto the central object. In the outer disk region there is a 
balance between inward advection (that wants to bring the outer footpoints to the central object), outward diffusion, and continuous field 
generation due to the action of the Cosmic Battery. The rate of growth of the accumulated magnetic field keeps decreasing and eventually 
saturate. 
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Figure 3. Evolution of the magnetic flux 4/ accumulated on the 
central object (in units of 4/ a cc) for various values of the magnetic 
Prandtl number (from top to bottom: Vm — 0.2, 0.4,1,10,100, and 
100 with finite optical depth in the disk). Time in units of t acc 

mulated magnetic field pressure B 2 / 47 r becomes compa¬ 
rable to the ram pressure of the flow). In Fig. 3, 

we have normalized to a characteristic value \l/ acc that 
corresponds to the amount of flux accumulated interior 
to risco generated by the Cosmic Battery electromotive 
force term in one accretion time, namely 


7rr i 2 sco i; CB(pisco)c 


Pr(nsco) 

W 0 \ Lcr T e 
v r 


Lcjt 


■, 10 14 G cm 2 
‘0.4 G Trrfgco 
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/ L \ ( M 


\L Edd J \10M & 


( 21 ) 


already at times t ~ 10 f acc when the inner and outer 
parts of the poloidal magnetic field loops begin to dis¬ 
connect (Fig. 2b,c). Beyond that time, the field lines of 
the return polarity which are anchored at larger distances 
on the disk are now wound by the disk rotation at their 
footpoints, and not by the differential winding between 
the inner and outer footpoints of the loop (which, as we 
said, are now disconnected). When this happens, the 
component of the return (outer) field line changes sign 
(becomes positive), as can be seen by the change of color 
near the surface of the disk. A dynamic magnetospheric 
current sheet forms between the two regions of opposite 
poloidal field polarity which in our very-low 77 case lies 
a few degrees above the surface of the disk (Figs. 2d-f). 
The ongoing (numerical) reconnection forms plasmoids 
that outflow fast along the current sheet. Plasmoid for¬ 
mation and reconnection in the current sheet may be 
responsible for particle acceleration and dissipation of 
energy. Magnetic flux is continuously advected on the 
central object, but at a continuously reduced rate (see 
Fig. 3). 

We have also run simulations in which we implement a 
long but finite optical depth in the disk according to the 
second term in the first line of eq. [5j for the same value 
of the disk diffusivity (low 77 ). This effectively kills the 
Cosmic Battery beyond a radial distance of a few nsco ■ 
In that case, the field saturation is much clearer as can 
be seen in Figs. 3 & 4. 

Notice that, in the very low 77 case, the saturation value 
of the magnetic flux that was brought to the central ob¬ 
ject is very low, thus, the observed current sheet activity 
observed in this case is probably astrophysically insignif¬ 
icant. The big difference from the standard case is that 
in the latter, the magnetic flux continues to grow at a 
steady rate, presumably up to equipartition. 


What is important to notice is that the magnetic field 
that threads the disk interior remains always at a very 
low value since, inside the ISCO it is continuously ad¬ 
vected toward the central boundary (the black hole hori¬ 
zon) , whereas outside it is continuously diffusing outward 
through the accretion disk. As long as the field remains 
at such small values, the MRI, the physical mechanism 
believed to be responsible for the disk turbulence (Balbus 
& Hawley 1992), does not saturate. 

4.2. The very low 77 case 

The field evolution in the very low 77 case is very dif¬ 
ferent. Beyond the accretion time t acc needed for matter 
to reach the central object, the mechanism begins to sat¬ 
urate. As we pointed out above (and as can be seen in 
Fig. 2a), at times up to a few times f acc the evolution is 
the same as in the standard case, namely the inner field 
footpoints are carried by matter falling to the central ob¬ 
ject and the outer footpoints are rotating with the disk. 
However, the accretion disk is no longer diffusive and the 
outer footpoints do not diffuse easily outward inside the 
disk. In fact, the accretion flow wants to bring also the 
outer footpoints to the central object. The whole loop 
would be swallowed by the central object, but this is not 
the case since new magnetic field loops are continuously 
generated by the Cosmic Battery in the optically thin 
part of the disk. Due to this, the mechanism instead 
saturates. We see the deviation from the standard case 


5. DISCUSSION AND CONCLUSIONS 

In the previous section we argued that the main fac¬ 
tor that allows for the steady growth of the accumulated 
magnetic flux is the value of the magnetic Prandtl num¬ 
ber Vm in the disk. We run several simulations for dif¬ 
ferent values of the magnetic Prandtl number in the ac¬ 
cretion flow (Fig. 3). In these simulations, other phys¬ 
ical parameters are kept constant. It is interesting to 
see that in all our runs the magnetic flux advected to 
the central object within a few accretion times t acc re¬ 
mains almost the same (on the order of 'P acc )• Within 
that short amount of time, the mechanism does not feel 
the diffusivity properties of the disk. As time goes by, 
only diffusive disks allow the magnetic flux to grow. As 
can be seen in Fig. 3, the rate of growth is closely cou¬ 
pled with the magnetic diffusivity of the disk. When 
the Prandtl number is low (high diffusivity), the rate of 
growth is bigger. For the lowest values that we have tried 
(V m = 0 . 2 ,0.4, 1 ), the growth of the magnetic flux is 
linear with similar slope after about f acc from the begin¬ 
ning of the run. As the Prandtl number increases beyond 
unity (which corresponds to lower and lower diffusivity), 
the rate of growth of the magnetic flux continuously de¬ 
creases with time, approaching zero asymptotically (no 
growth, i.e. saturation). There seems to be a critical 
value of the magnetic Prandtl number for this transition 
in the range 1 < V m crit < 10. For Prandtl numbers 
higher than 100 the situation is identical to the very low 
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Figure 4. Same as Fig. 2 but now an optical depth of order nsco is included in the disk and the Cosmic Battery dies out beyond a few 
times risco- Saturation is reached very quickly and beyond time t « 5t a cc no more magnetic flux is brought to the central object. 
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Pm = 0.4 
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Figure 5. Configuration (f) of Fig. 1 ("P m = 0.4, t = 38.4£ a cc) shown over a larger spatial scale. The central flux accumulation proceeds 
unimpeded. However, at larger distances, the disk develops complex dynamic magnetic activity. 


rj case. 

Our present realistic simulations of the Cosmic Bat¬ 
tery at the inner edge of the accretion disk around a 
black hole confirm that this is an important astrophysical 
mechanism that may indeed account for the generation of 
the large scale magnetic field that threads astrophysical 
jets and disks. Under standard accretion disk conditions 
(magnetic Prandtl numbers of order unity), the Cosmic 
Battery generates a large scale magnetic field on the or¬ 
der of B acc ~ 1 G inside the ISCO of a ten solar mass 
black hole within about one accretion time on the order 
of a few milliseconds. Assuming that the black hole is 
accreting at about 1% of its Eddington rate, in order for 
the accumulated magnetic field to reach an equipartition 
value on the order of 10' G, one has to wait for about 10 9 
accretion times, or equivalently for about t eq ~ 10 7 sec¬ 
onds (a few months). The above numbers obviously scale 
with mass as can be seen in Table 1. 

We have shown that the Cosmic Battery does not in¬ 
crease the value of the magnetic field that threads the 
accretion disk. The magnetic field remains several or¬ 
ders of magnitude below equipartition, and does not in¬ 
hibit the action of the MRI in the disk. Of course, as 
the magnitude of the accumulated field inside the ISCO 
approaches equipartition, the mechanism will saturate, 
and the accumulated field will probably escape outward 
through the disk due to magnetic Rayleigh-Taylor insta¬ 
bilities around its inner edge (Tchekhovskoy et al. 2010, 
Contopoulos & Papadopoulos in preparation). This is in¬ 
teresting because when that happens, one expects to also 


see magnetically driven outflows (winds) from the disk. 
In other words, unless the accumulated field approaches 
equipartition around the inner edge of the disk (and thus 
becomes unstable to magnetic Rayleigh-Taylor instabili¬ 
ties), we do not expect to see any strong wind emanating 
from the disk. This is contrary to what is expected in the 
scenario where the magnetic field is brought in from large 
distances through the disk, and thus threads the disk 
and generates strong steady disk winds. It is interesting 
that, although the overall growth of the generated and 
accumulated magnetic flux proceeds as described above, 
the return field that threads the accretion disk at large 
distance beyond nsco shows a very complex behavior 
that results from the interaction between the accretion 
flow, the force-free magnetosphere and the surface tran¬ 
sition region (Fig. 5). Such dynamic behavior is a far cry 
from the smooth steady-state configurations obtained in 
old calculations of the disk-wind interaction. Our results 
may have implications for the variability observed in flar¬ 
ing XRBs where a strong radio jet appears only when the 
system transitions from the hard to the soft high state 
(Kylafis et al. 2012). 

Similar complex variability is also observed in the nu¬ 
merical simulations of Parfrey et al. (2015) which are 
similar in concept to the ones presented in this work, 
although they solve only for the magnetospheric field. 
In the latter, the magnetic flux loops that inflate and 
open up are sustained by the disk turbulence as a bound¬ 
ary condition, and this is the reason highly variable re¬ 
connecting current sheets appear everywhere throughout 
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Table 1 

Timescales to reach equipartition 


M 

B a cc 

tacc 

Beq 

teq 

1 

0.1 

10 -3 

10 8 

one week 

10 

0.01 

10~ 2 

10 7 

a few months 

10 9 

10” 9 

10 ® 

10 3 

several Gyears 


Note. —M (in solar masses): the mass of 
the central black hole; B acc (in Gauss): the 
magnetic field accumulated inside the ISCO 
in one accretion time t acc (in seconds); B eq 
(in Gauss): the nominal equipartition mag¬ 
netic field; t eq : the time needed for the Cos¬ 
mic Battery to build B eq . 

their simulations. Our physical picture differs in that we 
do account for the origin of a large scale magnetic flux 
of a well defined polarity that threads the disk and the 
inner compact object. In our case, only one reconnecting 
current sheet develops at the interface between outgoing 
and returning magnetic field lines. 

Obviously, more investigation is required. Our work 
will be extended soon with general relativistic simula¬ 
tions that take into account the full distorted radia¬ 
tion field generated by the innermost hot accretion flow 
around the central spinning black hole (Koutsantoniou 
& Contopoulos, in preparation). 

This work was supported by the General Secretariat 
for Research and Technology of Greece and the European 
Social Fund in the framework of Action ‘Excellence’. 
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